Stochastic potential and quantum decoherence of heavy quarkonium in the 

quark-gluon plasma 



Yukinao Akamatsu^ and Alexander Rothkopf^ 

Kobayashi-Maskawa Institute for the Origin of Particles and the Universe (KMI), Nagoya University, Nagoya 464-8602, Japan 
'^Fakultat fur Physik, Universitat Bielefeld, D-33615 Bielefeld, Germany 
(Dated: March 23, 2012) 

We propose an open quantum systems approach to the physics of heavy quarkonia in a thermal 
medium, based on stochastic quantum evolution. This description emphasizes the importance of 
collisions with the environment and focuses on the concept of spatial decoherence of the heavy 
quarkonium wave function. It is shown how to determine the parameters of the dynamical evolution, 
i.e. the real potential and the noise strength, from a comparison with quantities to be obtained from 
lattice QCD. Furthermore the imaginary part of the lattice QCD heavy quark potential is found to 
be naturally related to the strength of the noise correlations. We discuss the time evolution of QQ 
analytically in a limiting scenario for the spatial decoherence and provide a qualitative 1-dimensional 
numerical simulation of the real-time dynamics. 

PACS numbers: 



I. INTRODUCTION 

The fate of heavy quarkonium states {QQ) at very high 
temperature is a long standing puzzle ever since suppres- 
sion of J/vt (the cc ground state with = 1 ) has been 
proposed as a prime signal for the formation of the de- 
confined state of QCD matter, the quark-gluon plasma 
(QGP) [HE]. The production of J/* and T {hh states 
with = 1~) in heavy-ion collisions has been measured 
in detail at the Relativistic Heavy Ion Collider (RHIC) [5] 
and more recently at the Large Hadron Collider (LHC) 
[HH]. Partial suppression of J/^* j3j and the relative sup- 
pression of excited states such as T(2S, 3S) in comparison 
to T(1S) [5] have indeed been observed. 

In order to attack the question of QQ survival, sev- 
eral approaches have been deployed. Maximum Entropy 
(MEM) analyses [6] of J/^P and T spectral functions in 
lattice QCD have revealed a possible survival of these cc 
and hh bound states up to around T ~ 2Tc with Tc being 
the critical temperature for the hadron-quark transition 
[7j. The observed spectral structures are interpreted as 
a sign of J/'^ survival by some authors T], while for 
others they represent mere threshold enhancement 
The unclear definition of what constitutes a heavy quark 
bound state is clearly a drawback of this otherwise solid 
approach. 

On the other hand, purely real potential models, based 
on quantities, such as the difference in the free or inter- 
nal energies [2| of a medium with and without a QQ in- 
serted, are struggling with the absence of a Schrodinger 
equation derived from QCD. This leads to ambiguities 
in their definition of a heavy quark potential. Progress 
has been made in deriving a Schrodinger equation from 
QCD by evaluating the late time evolution of the real- 
time thermal Wilson loop [TUlI12j . An initial perturba- 
tive study at very high temperatures, based on the Hard- 
Thermal-Loop approximation, yielded a complex poten- 
tial, whose imaginary part is induced by collisions with 
the light plasma particles. Recently [TJ] showed how to 



relate the potential to the spectral decomposition of the 
thermal Wilson loop and based on lattice QCD simula- 
tions confirmed the existence of an imaginary part above 
Tc. 

In this paper, we instead propose an open quantum 
systems approach (see [111 [15] for earlier attempts), using 
a stochastic, i.e. a fully dynamical description of QQ in 
the QGP. The noise inherent in our stochastic treatment 
is a result of interactions between the QQ open system 
and the medium, which is traced out from the descrip- 
tion. It not only naturally gives a physical meaning to 
the imaginary part found in the heavy QQ potential from 
lattice studies but also provides an indispensable ingre- 
dient in discussing the time evolution of QQ, i.e. spatial 
decoherence of the wave function due to the collisions 
with plasma particles. 

Without spatial decoherence, such as in real poten- 
tial models, the effects of collisions with the environment 
are underestimated, whereas having an explicit imagi- 
nary part in the potential will lead to an unabated sup- 
pression of all available states. Both of these situation 
are expected to be unphysical and we will show how our 
proposed model fills this gap in the following sections. 

Our investigation makes close contact with recent lat- 
tice QCD studies [T3| by relating the strength of the noise 
to the spectral function of the thermal Wilson loop, thus 
making possible the determination of the required pa- 
rameters at all temperatures non-perturbatively. 

This paper is organized as follows. In section[ll| we for- 
mulate the stochastic equation of motion for the heavy 
quarkonium wave function and give a physical interpre- 
tation of the correlation present in the noise term. The 
limitations of our proposed formulation are listed and dis- 
cussed. In order to connect the stochastic model to ac- 
tual physics, we show in Sec. |III| how to determine some 
of the basic parameters of the stochastic dynamics from 
first principles lattice QCD simulations. In Sec. |IV[ we 
derive a master equation for the density matrix of heavy 
quarkonium states from the underlying stochastic time 
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evolution. We then trace out the global motion from the 
density matrix, which enables us to reduce the two-body 
problem to a one-body problem. In Sec. |Vj we first hst 
several criteria for bound state survival and determine 
their values in one dimensional numerical simulations of 
the stochastic dynamics based on different parameter sets 
adapted from the literature. Finally in Sec. |VI[ we sum- 
marize our work and give an outlook toward how to arrive 
at a more complete understanding of the heavy quarko- 
niimr states in the QGP. 



II. STOCHASTIC EVOLUTION OF AN OPEN 
QUANTUM SYSTEM 

Suppose we insert a heavy quark bound state into the 
thermal medium and let it evolve in time. In a situation 
where a classical description is adequate, the relevant 
dynamics of the heavy particle is governed by Langevin's 
theory of Brownian motion. On the other hand, in a sit- 
uation where quantum effects play a crucial role, the dy- 
namics should be described by the theory of open quan- 
tum systems [TQ. In the latter case, interactions with 
the light particles of the heat bath lead to fluctuations in 
the time evolution of the heavy particle wave function, 
or in other words, to fluctuations of the Hamiltonian, i.e. 
the potential. In the phenomenologically relevant region 
of the strongly coupled quark-giuon plasma just above 
the deconfinement phase transition, surely the latter ap- 
proach is required. 

The full quantum system is described by a density ma- 
trix and its time evolution is governed by the full Hamil- 
tonian (in the Schrodinger picture): 



function ^qq{X, t). Time evolution is thus described by 
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with some initial condition p(0) — Psys® Pmcd- In the fol- 
lowing, we are only interested in the physics of the sub- 
system, which is described by a reduced density matrix 
Psys{t) = Tvmed{p{t)}- The dynamics of Psys{t) is often 
unraveled into a stochastic evolution in Hilbert space. It 
can be described e.g. by a master equation of Lindblad- 
form 17J, corresponding to stochastic quantum evolution 
that incorporates quantum state diffusion |18j or quan- 
tum jump processes [19]. 

In our approach we first wish to formulate a wave func- 
tion based description of the stochastic time evolution of 
a non-relativistic heavy quarkonium state, whose pair an- 
nihilation is neglected. We use the short hand notation 
X = {xi,X2) to denote the position of the constituent 
quarks, written as single vector in a 6-dimensional space. 
The corresponding distance in this vector space reads 
AX = \X - X'\ ^ y/{xi^ x[y + {X2- x'^)^. The 

(X) 

stochastic time evolution operator Uq {t\0) has to be 
unitary and we assume that the dynamics it describes is 
both Markovian and linear in the heavy quarkonium wave 



^QQix,t) = ui''\tmQQix,o), 
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where T denotes the time-ordered product and the 
stochastic term in the Hamiltonian Q{X,t) corresponds 
to Gaussian white noise with the characteristics: 

{e{x.t))^o, {eix,t)eix',t')) = hTix,x')Stt'/At.{A) 

Here we choose as a time scale At, during which the 
plasma particles experience many collisions but the heavy 
quarkonium state does not change considerably. Such in- 
termediate time scale At is expected to exist due to the 
hierarchy in the mass scale M ^ T. The fluctuation of 
the Hamiltonian derives from the variations in the trans- 
ferred energy during At originating from collisions with 
the medium particles. The squared strength of the fluctu- 
ations at X, which is given by hr{X , X) / At, represents 
how frequently and effectively the energy transfer takes 
place between the heavy quarkonium and the surround- 
ing medium. In a uniform and isotropic system, which we 
consider, T{X,X) depends only on \xi — X2\- Since the 
relativistic medium particles are of size Ith = iTrhc/kBT, 
the fluctuations, or equivalently the transferred energy, 
is expected to be correlated at X and X' if AX < /th- 
Therefore we expect that T{X,X') '--^ if AX is larger 
than ~ /th- 

We take as cutoff length scale for this description the 
Compton wavelength of the heavy quarks /cut = h/Mc <C 
/thj so that the non-relativistic treatment is valid. Since 
the typical momentum of heavy quarks near thermal 
equilibrium is y/ksMT = 27rS/ZQ_th, there is a hierarchy 
lent ^ ^Q.th ^ 'th- It is this hierarchy that requires the 
effective stochastic dynamics of heavy quarks to be uni- 
tary, because collision processes have almost no chance 
to excite heavy quarks to a very high momentum state 
beyond the cutoff momentum scale 27r/i//cut-Note that 
this argument implies that a mechanism for thermaliza- 
tion must also be incorporated into a complete effective 
dynamics of heavy quarks. 

Given the manifestly unitary evolution in Eq.(|3|, the 
stochastic differential equation for the wave function can 
be written down by expanding Uq {t+At\t) in At, yield- 
ing up to 0(Ai3/2) 



Ui''\t + At\t) 

« 1 - ^ {H{x) + eix,t)} - ^e(x,t)2 



(5) 
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Here we define a new complex noise 5(X,t) in terms of 
e(X, t) so that t)) = 0. The stochastic differential 

equation for the wave function in Ito discretization thus 
reads 



(6) 



= [h{x) - '-r{x,x) + E{x,t)}^QQ{x,t). 

Note that the expansion of the unitary time evolution 
operator has induced what appears to be a non-unitarity 
in the Schrodinger equation. Since the origin of this con- 
tribution however are the noise terms of Eq. ([5| , it is only 
after taking the ensemble average that damping of the 
wave function is observed. 

We would like to remark that the formulation of the 
stochastic dynamics based on Eq.([3| is most likely not 
able to describe the quantum analog to the drag force in 
classical Langevin theory. Classically, this corresponds 
to an absence of a friction term for momentum diffusion, 
i.e. the momentum distribution of the heavy quarks dif- 
fuses without resistance. The effects related to these ab- 
sent terms become important to the dynamics once the 
typical energy of the heavy quarks reaches and exceeds 
the temperature. Hence, roughly only for time scales 
smaller than the heavy quark relaxation time our descrip- 
tion should be valid. Note that heavy quark relaxation 
time can be much longer than that of the medium At by 
a factor ~ M/T, so that the applicability of our descrip- 
tion is not as severely limited as it might appear at first 
sight. 

The above derivation of Eq. ^ is based on the expan- 
sion of Uq {t + At\t) in At. However, since the time 
step is related to physical scales as At ~ h/g'^ksT, it 
may not be an expansion in small dimensionless number. 
This subtlety is overcome by introducing N independent 
Gaussian noise variables Qi{X,t) (i — 1, • • • ,N) which 
satisfies Eq.Q with At replaced by At/N. If we define 
Q{X,t) = j^T.f-^i{X,t), Q{X,t) satisfies the origi- 
nal definition Eq.Q. The new noise Oi(X,i) has a nat- 
ural interpretation as fluctuation of the Hamiltonian in 
Eq.Q with artificial discretization time scale At/N . We 
then assume that Eq.([3| can still be expanded as time- 
ordered product of white noise, now with respect to the 
Oi(X,t)'s. Taking the noise to be independent even on 
time scales smaller than At is a prerequisite to obtaining 
a Markovian master equation of the form present in Sec- 
tion IV. As we observe heavy quarkonium on time scales 
longer than At, we discretize Eq.([3]), expand Eq.([5]), and 
simulate Eq.Q on a smaller timescale At/N instead of 
At. In the remaining part of this paper, we distinguish 
dt = At/N from At. 



III. SPECTRAL ANALYSIS 

Now that we have formulated the stochastic evolution, 
one has to ask how the unknown functions V{X) and 



T{X,X') in Eq.lll can be determined. Since we are 
interested in describing physics in the strongly coupled 
region of the quark-gluon plasma, we are urged to find 
non-perturbative means to do so. 

In analogy to the concept of the Nambu-Bethe- 
Salpeter (NBS) wave function at zero temperature, the 
averaged wave function qq[X ,t) of heavy quarkonium 
at finite temperature is identified by a suitable matching 
prescription with the following gauge invariant mesonic 
correlator (t > 0): 

D>{X,t) = (M(X,t)Mt(Xo,0))T = (*QQ(X,t)),, (7) 

where M'^^ {X, t) is the point-split heavy quarkonium an- 
nihilation (creation) operator and {■)t denotes the field 
theoretical average over the medium degrees of freedom 
at finite temperature [101 US]- In the case of infinitely 
heavy quark masses, where mq ^ T and mg 3> Aqcd, 
the left hand side of Eq. ([7| can be calculated in terms of 
the thermal Wilson loop 



D>^^^{X,t)^{W{X,t))T^{e^V 



(8) 



where A^{x) denotes the gauge field of the medium glu- 
ons, which is integrated along the rectangular contour 
connecting Xq at t = and X at t. The right hand side 
of Eq.([7]) becomes 



(*QQ(^,i)>H «exp 



-fW)4r(x,x) 



(9) 



In order to continue we have to remember that all equa- 
tions appearing so far are expressed in real-time t. As 
we wish to use Monte-Carlo simulations of lattice QCD, 
formulated in Euclidean time to evaluate the Wilson loop 
non-pcrturbatively [20.^, their results need to be analyti- 
cally continued. One possible strategy is to use spectral 
functions, defined from the forward and backward prop- 
agator 



a{X,uj) ^ D>{X,uj) - D<{X,uj) 



(10) 



Since the backward correlator for infinitely heavy quarks 
vanishes D^^^^{X ,t) = 0, one finds that the spectral 
function a{X,ijj) is nothing but the Fourier transform of 
the forward correlator D-^{X, t) and thus of the real-time 
Wilson Loop. Analytic continuation yields 



W{X,t) 



doj e-"^CT(X,a;), 



(11) 



which tells us that the relevant spectral information can 
be extracted from imaginary time data accessible on the 
lattice. In practice the numerical determination of the 
above spectral function from the Euclidean Wilson loop 
is based on a form of Bayesian inference called Maximum 
Entropy Method (MEM) [6j. 
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In case that the spectrum of the Wilson loop exhibits 
well defined structures, e.g. peaks of Breit-Wigner form 
[H], the real potential V{X) and the local noise corre- 
lations r(X, X) of Eq.(|9]) are readily obtained from the 
position and width of the lowest lying peak [T^ : 



a{X, Lu) cx 



r{x,x)/2 



{huj~v{x)y + {r{x,x)/2) 



2 • 
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In the context of our model, the width of the spectral 
function is interpreted as uncertainty in the actual value 
of the real potential in the open system. Consequently 
time evolution is always unitary for each realization of 
the stochastic ensemble. In this sense it characterizes the 
decoherence between the initial and current wave func- 
tion. This is in contrast to the deterministic notion of 
the width as an imaginary part, i.e. an explicit contribu- 
tion to the Hamiltonian, which violates hcrmiticity and 
leaves open questions about conservation of the proba- 
bility density. 

Even though we have already obtained two main in- 
gredients of our model, the real potential and the local 
noise, we are still missing another piece of the puzzle, 
i.e. information on the off-diagonal noise terms. A quan- 
tity that might allow extracting this set of parameters is 
introduced in the next section. As a first ansatz, to al- 
low simulations based on Eq.([6]), one might assume that 
such spatial correlations within the fluctuations, i.e. the 
AAT-dependence of T{X,X'), will decay with the ther- 
mal wavelength as exp{—{AX)'^/l^^). 



IV. THE MASTER EQUATION 

In principle, one can simulate the stochastic time evo- 
lution in Eq.([6]) directly, but in order to precisely define 
the time evolution of the heavy quark bound state in a 
thermal medium, we are urged to go one step further in 
our analysis. To this end we derive a master equation 
for the reduced density matrix of states pqq{X, X',t) = 
{^QQ{X,t)^*Qq{X' ,t))-. Noting that an infinitesimal 
time step in the time evolution of Pqq{X, X' , t) is given 
by 



PQQi^;X',t + dt) 



(13) 



= {Ui''\t + dt\t)U^^ >*{t + dt\t))^ PQQ{X,X',t), 

and calculating {U^\t + dt\t)U^'^*{t + dt\t))- using 
Eqs.Q and ([S]), it is straightforward to derive the fol- 
lowing relation: 



AX') 



dpQQ{X,X',t) 
dt 



H{X)- H{X') 
ih 

, F{X,X') 



PQQ{X,X',t) (14) 

X,X',t), 



F{X,X') = T{X,X') 



T{X,X) + TiX',X') 



The conservation of the trace Tvx{pqq{X , X' ,t)} = 
J d^XpQQ{X,X,t) — 1 can be proved easily. Since the 
time evolution of this averaged quantity now depends also 
explicitly on the off-diagonal spatial correlations in the 
noise, it appears as a promising candidate to extract the 
yet undetermined values of T(X' , -^')- Future work thus 
needs to focus on expressing the density matrix of states 
for the heavy quark system as field theoretical operator 
amenable to lattice QCD simulations. 

An equivalent master equation is known from the stud- 
ies of scattering models and applications of influence- 
functional techniques to random potentials [22) . Its ap- 
plicability to the physics of heavy quarks, however, has 
so far not been exploited. Note that in comparison, the 
recently proposed quantum optical master equation |15| 
is applicable where the rotating wave approximation is 
valid, i.e. in discussing transitions between deeply bound 
states [IB] . Our approach on the other hand is geared to- 
ward the description of loosely bound states and melting 
phenomena, since it explicitly describes both constituting 
particles individually. 

As we are interested in the relative motion of the 
heavy quarks we trace out the global motion. Defin- 
ing the further reduced density matrix PQQ{r,r' ,t) = 

Trfl{pQQ(^«,r,^i?',r',i)} = ! 'PRpQQ{^R,r,XR^r' -.t), 

where = {R+r /2, R—r /2), we obtain from Eq. ( 14 ) 
the master equation for pQQ{r,r',t): 



dpQQ{ry,t) _ h{r)-h{r') 



ih 

f{ry) 



PQQir,r',t) 
-j{r,r\t), (15) 



dt 



where h{r) = 2M — ^ ^['' +v{r) and f{r,r') = ^{r,r') — 
{j{r,r) + j{r' ,r')}/2. Here we define v{r) = V{XQ,r) 
and 7(r,r') = r(Xo,r, -X^o,r') using translational and ro- 
tational invariance. Correspondingly, the stochastic dy- 
namics of the relative coordinates of the heavy quark pair 
is given in Ito discretization by 



d 



(16) 



[h(r) - ^7(r, r) + ^(r, i)}^QQ(r, t), 



idt 



with {e{r, t)) = and {e{r, t)9{r', t')) = h-y (r, r') Sw/dt. 

In order to discuss the time evolution of heavy quarko- 
nium in the next section, we can take two equivalent, 
while differently nuanced standpoints. From the theoret- 
ical side, one wishes to describe the time evolution of the 
system by as simple equations as possible, which is best 
done in the basis of eigenstates of the Hamiltonian gov- 
erning the dynamics h(r). On the other hand to build a 
bridge to experiment and its notion of QQ suppression, 
one asks how much of an initial bound state has survived 
after a certain time t. This question is best answered 
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when the evolution of the system is observed in the ba- 
sis of the eigenstates of the initial vacuum Hamiltonian 
/i™'^(r) that governs the bound state. 

We thus transform Eq.(|15|) by expanding 



PQQ{ry,t) =J2cZAt)Mr)r„^{r') (17) 

nm 

= T.<mimn{rWmir') (18) 

nm 

in a complete set of eigenfunctions ipnif) (with eigenvalue 
e„) of h{r) or (j)n{r) (with eigenvalue e„) of /i™'^(r). The 
expansion coefhcients are thus given by 



(19) 

In terms of the in- medium state admixtures c™„, the time 
evolution of the system according to Eq.( 14 ) takes on the 
following form: 



(t ) ^ — —c 



in n 



kl 



]^Y.{^r.kctAt)+C^:,(t)^Ura}, (20) 



2h 



where 7„fcjTO EE / d^rd^r'^p:^{r)^p^ {r')j {r,r') ipkMiJmir') 
and 7„™ = / d^rtp;^{r)'y{r,r)tpm{r). Here 
c™^(t) and 7„m are hermitian matrices and 

lnk,lm = llm,nk = "/kn,ml- 



V. STOCHASTIC EVOLUTION OF HEAVY 
QUARKONIUM IN THE QGP 

To model the physics of heavy quarkonia related to rel- 
ativistic heavy ion collisions at RHIC and the LHC, we 
use the following setup. The initial heavy quarkonium 
is taken as the ground state of the vacuum Hamiltonian 
^vac ^^•■j ^ ciiaracterized by a Coulombic and linearly rising 
potential at small and large separation distances respec- 
tively. Once such a bound state enters the region in which 
a quark- gluon plasma is present, the interactions with 
this hot medium are taken into account through time 
evolution as given by Eq.(16) with appropriate model 
parameters v{r) and 7(r, r') corresponding to tempera- 
tures T > Tc ~170MeV. Depending on whether we wish 
to focus on the time evolution of the in-medium states or 
the initial bound-states in the QGP, we may use either 
the medium or vacuum bound state survival probabilities 



l<n<Af;" 



(21) 
(22) 



where TV™ (^b)is the number of bound states of h{x) 
(/i™'^(a;)). 



A. Analytic Considerations 



Among three different representations of the stochas- 



tic dynamics in Eqs.( 15 1,( 16 ), and (201, the qualitative 



feature of P™(f) can be discussed most clearly by using 
Eq.(20). The following arguments hold also for P^{t) as 
long as h{r) = ft,™''(r). 

In the simplest case where j{r,r') — 7, we can show 
that 7„fc,/„i = jSnkSim, jnk = jSnk, and thus c™m(0 = 
j(i), P™(i) = 0. Therefore with spatially uni- 



form fluctuations, quarkonium survival is completely de- 
termined by h{x), which is the case for the purely real 
potential models. 

On the other hand, in a case where the spatial correla- 
tions factor according to 7(r,r') = 7exp(— (Ar)^/2Z^Qj.j.) 
with Ar = |r — r'l and Zcon- ^ ^cut, we can show that 
Ink = J^nk and 7„fc,/„i as /corr/'cut -> 0. The latter 
relation is obtained by using the fact that the eigenfunc- 
tion ipnix) varies over length larger than Icut'- 



InKim ~ 1 ] d'rd'r'rn{r)ri{r)e Mr)i^m{r) 

< 7(v^^corr/Ut)^ ^ (U-r//cut ^ 0).(23) 



Then we obtain c™„(i) 



At) 



7 „m 
h nm 



l<n<NZ 



P™{t) = — ^P™(i). Therefore with very localized fluctu- 
ation, the damping is determined solely by 7 = 7(r,r). 
Since the correlation length ^corr of the fluctuation is be- 
low the cutoff scale Icnti such fluctuations excite wave 
modes which are beyond the model description, i.e. 
all the diagonal elements of c™„(i) decay as ^ e"^*/'' 
and Trr{/o(r, r', /:)} is not conserved. In the physical 
case ^corr = ^th; the above arguments suggest that the 
states with typical scale larger than ^th decay while those 
smaller than Zth remain undisturbed, supporting the ar- 
gument that the fate of bound states in essence is a 
hadronic thermometer. 



B. Id Numerical Simulations 

For a first qualitative exploration of the above ideas, we 
carry out one-dimensional numerical simulations of the 
stochastic dynamics at T ~ 2.33rc, based on different 
parameter sets adapted from a lattice QCD study [T3] 
(models (A) and (C)), the perturbative (PT) complex 
potential calculated in [10] (model (B)) and the classic 
Debye screening scenario of [T] (model (D)). (Throughout 
this section, we adopt the natural unit h = c = ks = 
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dx 
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dx 
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To accommodate the full heavy quarkonium wave 
function Eq.(16) is discretized on a line of iV = 512 
points, that covers the physical distance between x G 
[— 2.56fm, 2.56fm]. The initial wave function is deter- 
mined as eigenstate to the Cornell potential = 
— 1^ + ij\x\, with vacuum parameters a = 0.1 and 

cr = (0.4GeV)^. To implement the effects of dynamical 
quarks, i.e. to account for string breaking, we flatten the 
linear rise smoothly to a constant value at Tsh = 1.5fm. 
This choice leads to a total number of A''^ = 4 bound 
states. All survival probabilities P^{t) shown in this sec- 
tion are calculated based on the admixture of these four 
vacuum states contained in the system wave function. 

The stochastic dynamics are implemented using a 
Crank-Nicholson scheme with time step dt = O.OOOlfm. 
There the mass of the constituent quarks is set to M — 
l.lSGeV in order to be larger than the temperature, 
which we choose to be T = 2.?,iTc - 0.4GeV. In the 
first three runs (Figj2] and Figjs] left) , off-diagonal noise 
terms are set to zero, leading to an artificial correlation 
length Zcorr ~ dx < Zth = 15GeV^^. This yields a noise, 
which is stronger than expected from the physical value 
of the temperature. 

In the third example (Figjsjleft) we approximate the ef- 
fects of nonvanishing off-diagonal correlations 7(0;, a:') oc 
exp(— |a; — x'p/^^Qjj) by initializing noise in Fourier space 



{e{p, t)) = 0, t)e{p\ t')) = k{p)5 {p + p') Stt'/dt 

K{p)= J dxe^P-exp(-xVzLr) (24) 

with a normalized Kernel k{p) at each time step I23J. 
The strength of the diagonal correla tions is then set by 
simple multiplication with \/"f{x, x) after transforming 
back to coordinate space. 
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FIG. 1: Norm and Energy for the parameter set of model 
(B). We find unitary evolution, so that the wave function 
norm I'^qqC^)! = 1 and its energy remains purely real. After 
averaging, decoherence induces a dampening in \{tpQQ}e{t)\ 
as expected from Eq.([9]|. The diflterence in E{t) and E^{t) is 
dominated by the noise contribution to the system energy. 



Before we start a comparison of the different results 
according to the parameter sets (A)-(D), let us briefly 
look at generic features of the simulation from the view 
point of wave function norm 



IV'qqWP^ / dxil;*r,{x,t)tljnQ{x,t) 



(25) 



and energy 



E{t) = J dxi)QQ{x,t)h{x)ipQQ{x,t) 
E^t)^ j dxrQQ{x,t)K"^-{x)i,QQ{x,t) (26) 

Both merits and the limitations of the stochastic model 
of Eq.Q are already visible in these quantities, which we 
display in Figjl] based on the parameters of model (B). 

One finds that indeed unitarity is preserved for each 
individual run of the stochastic ensemble since the wave 
function norm remains at unity. On the other hand, per- 
forming the ensemble average according to Eq.([9]) intro- 
duces a dampening of K^'qq)!, which is readily observed 
in the decreasing dashed line in Figjl] 

Since our stochastic quantum evolution does not in- 
clude the quantum counterpart of the drag force, mo- 
mentum of the heavy quarks diffuses in momentum space 
without resistance by a friction term. Consequently one 
expects that the overall energy of the system over time 
will increase artificially and indeed one observes a linear 
rise in the upper two dashed lines in Fig{T] at intermedi- 
ate times. The artificially strong noise pushes the values 
of the energy to quickly exceed the non-relativistic cut- 
off ^ Af ^ l.lSGeV and saturates at late times at 
E ^ (7r/dx)2/M SOOOGeV due to the momentum cut- 
off on the lattice 7r/dx. To eliminate the unphysical linear 
rise in energy and to allow the system to thermalize even- 
tually will require a modification of our approach that 
takes into account the correct momentum dissipation, i.e. 
friction as e.g discussed in [21]. Although we do not know 
how the effect of friction changes the time evolution of 
the individual heavy quarkonium wave function in detail, 
we expect these neglected effects to become important 
once the energy reaches the temperature of the medium 
(i.e. E{t) r^T 0.4GeV). Thus the late time behavior 
observed from the stochastic dynamics presented in this 
study is to be understood within this limitation. 

Note that especially the simulations of model (A) and 
(B) are performed with rather short correlation length 
(27r/^corr ^ 125GeV). Therefore although qualitatively 
similar behavior should be observed in simulations with 
longer correlation length (Zcorr 27r/T), the rapid evolu- 
tion presented hereafter must not be understood as rep- 
resenting a realistic evolution time scale. 

Model (A): The first run shown on the left of Fig. 
1 is based on a parameter set that models thermal ef- 
fects as fully incorporated into the noise strength. The 
real potential v{x) = u™'^(a;) in this case is not modified 
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FIG. 2: One-dimensional model: (top left) Vacuum admixtures cj^„(t) of heavy quarkonium bound states and their survival 
probability P^{t) from the stochastic model, based on a lattice QCD inspired parameter set [13]. In the bottom left panel we 
show the corresponding values of the initial potential u™'^(a;), the real potential v{x) governing the dynamics and the diagonal 
noise strength 7(0;, a;), (energy levels of the initial bound states et as well as for resonances and continuum solutions e^.c to 
v^^'^{x) are given as reference once). Decoherence induces an almost exponential decay of the initial heavy quarkonium ground 
state at short times before a stable plateau is reached. We find that all available eigenstates (pn of v^'^'^{x) are excited and at 
late times their relative abundances all coincide, (right) Same quantities with a parameter set adapted from a perturbative 
(PT) study [To] where at T = 2.33Tc a Debye screened real potential v{x) = Re[i'^^(a;)] is accompanied by a small but finite 
noise term 7(2;, a;) = 2lm[v^'^ (x)]. 



from its T = form, while the strength of the fluctua- 
tions exhibits a Unear rise equal to that of the real part 
7(x, x) — a\x\. 

One finds that the medium via the noise, exponentially 
decreases the probability to find the initial ground state 
at early times t < 4fm. It is the concurrent excitation 
of higher lying states and their feed-down that at later 
times t > 7fm leads to a stabilization of the ground state 
admixture and thus halting its decline. 

The artificially small correlation length of the noise 
also appear to render insignificant the energy differences 
between the individual bound states. Hence already at 
times t > 2fm the noise induced admixtures of higher 
lying states become of equal magnitude. 

Following the time evolution to late times t > lOfm 
a constant survival probability P^{t — 00) w 10^^ 
emerges. We however refrain from attaching physical 
meaning to this feature, since the inclusion of the fric- 
tion can change the outcome at these times significantly. 

Model (B): A qualitatively different setup is shown 
on the right of Fig.l. Its parameters are obtained from 
an evaluation of the heavy quark potential v^'^{x) in re- 
summed perturbation theory (TU] . Here the effects of the 
medium interaction are present both as modification of 
the real potential towards a Debye-screened form 

-mo\x\ 

v{x) = Re[v^'^{x)] cx -niD n — (^7) 



with Debye mass mjj = 0.92GeV, as well as the presence 
of a noise term 

, X „T r PT/ M f°° dz z sin[zmn.Tl \ 
7(a:,x) = 2Im[t;PT(a,)]« / -^—-Al ^^^V 

Jo (2 + Ij"^ V zmnx / 

(28) 

In comparison to model (A), the strength of the noise 
7(x, x) as well as the real potential are much weaker. We 
find that in total it takes longer for the ground state to 
become suppressed. Note that in the absence of noise, 
the parity even initial ground state is able to mix only 
with other parity even states under the time evolution 
of the P-even Hamiltonian h(x). The relatively strong 
population of the state 03 compared to the lighter (^2 
and (/14 thus tells us that it is such mixing processes that 
are dominant at times t < 5fm. 

Only after this point in time do the decoherence in- 
duced parity odd states contribute to similar strength as 
(1)3. At late times the asymptotic value for the suppres- 
sion P'^{t — > 00) w 10~^ in model (B) is very similar 
to the value observed in model (A), even though its real 
potential as well as the diagonal noise are very different. 

Model (C): Let us thus take a look at how spatial cor- 
relations in the noise influence the suppression pattern. 
Purely local noise corresponds to a very high medium 
temperature if interpreted as a correlation length of the 
size of the lattice spacing. Thus a realistic description 
around the deconfinement phase transition will need to 
allow for larger values of Zcon- ~ ^ ■ On the left of Figp5^ 
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FIG. 3: (left) Same u(a::) and ^(x,x) as in model (A) but with Gaussian damped spatial correlations 7(3;:, a;') 
e noif 

ao|x| 



exp[— |x — 

x'l'^/lcorr] induced in the noise terms using /corr ~ 4dx. (right) The classic Debye screening scenario based on a completely real 



potential vix) 



with the Debye mass chosen as mB = 5GeV. 



we show the results for a noise with the same local linear 
rise as in Figj2] left, but with Gaussian damped spatial 
correlations 7(2;, a;') cx exp[— |x — x'^/l^^j.^.] induced by a 
convolution with an appropriately normalized Kernel in 
Fourier space using /con- = 4dx. 

One finds that even though the integrated strength of 
the noise is the same, i.e. ^ k(p) — 1, its effects are 
not as localized anymore and thus the suppression of the 
ground state proceeds slower. The pattern appears to be 
different also in a qualitative way. We observe e.g. that 
the excitation of the state ^3 does not anymore show the 
overshoot present in all the other scenarios above. 

With a noise corresponding to a lower apparent tem- 
perature, the difference in mass, i.e. eigenenergy, be- 
tween the individual bound states becomes more impor- 
tant. Fig|3] (left) shows at small times t < lOfm that the 
amount of population of the states appears to be ordered 
according to their masses. From the higher lying states, 
the lightest is consistently excited most strongly. It is 
the relative abundances of the heavy quarkonia states 
that are closely related to the off-diagonal components of 
the noise correlations. 

Model (D): Taking the route from Model (A) via 
model (B), we arrive at the Debye screening scenario [I], 
where all medium effects are modelled as a modification, 
i.e. a weakening, of the real potential. If real and imag- 
inary part were to each encode part of the thermal ef- 
fects on the QQ system, the medium effects that were 
described in model (B) by noise should now be captured 
by an increased amount of screening present. Hence we 
choose tud = 5GeV. 

As shown on the right of Fig|3] one finds that indeed 
mixing of the eigenstates of h'^^'^{x) occurs and the initial 
bound state becomes suppressed. P^{t) however does not 



decrease exponentially not even monotonically. 

Note that populating states other than the initial 
ground state in this case can only be facilitated by mix- 
ing through the Hamiltonian. This leads to the observed 
pattern of the cafe's, where only parity even states are 
present. (Note that we plot all P-even states up to <j)s) 

So far we can summarize our findings: 

• The presence of noise leads to spatial decoherence 
and contributes at short times to an exponential 
suppression of the initially present heavy quarko- 
nium ground state. 

• Decoherence also populates the higher lying states. 
Its effect can differ significantly from the mixing in- 
duced by the Hamiltonian h{x) if it is not restricted 
by selection rules such as e.g. parity. Feed-down 
from these states leads to a state with respect to 
P^(t) at late times. 

• Details of the time evolution, such as suppression 
speed and the population ratios of heavy quarkonia 
states already at early times appear to be directly 
related to the structure of the off-diagonal noise 
correlations. 

• The overall suppression at late times in terms of 
P^(t) is quite insensitive to the details of v{x) and 
7(2;, x'). This behavior is however an artifact due 
to neglecting the effect of friction, which becomes 
important at late times. 
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VI. SUMMARY AND OUTLOOK 

In this study, we have proposed an open quantum sys- 
tems approach for the description of QQ states in a ther- 
mal medium using stochastic evolution. The merits of 
the stochastic evolution are the possibility 

• to give a dynamical description of the non- 
relativistic heavy QQ evolution in the quark gluon 
plasma 

• to give a physical meaning to the existence of an 
imaginary part in previous studies of the heavy 
quark potential 

• to make accessible the concept of spatial deco- 
herence in a potential based description of heavy 
quarkonium suppression 

while its applicability is limited to time scales shorter 
than the heavy quark relaxation time due to 

• the absence of the friction and the resulting inabil- 
ity to thermalize, expressed in a linearly rising en- 
ergy lim^^oo^W ^0. 

We have shown how to relate the basic model parame- 
ters governing the dynamics, such as v{r) and 7(r, r), to 
the spectral decomposition of the thermal Wilson loop, 
which in turn can be determined from lattice QCD non- 
perturbatively. It is found that the presence of noise 
leads to spatial decoherence and the consequent damp- 
ening of the QQ wave function, which plays a central 
role in describing the time evolution of the heavy quark 
bound states. One-dimensional model calculations based 



equation, as we have done in deriving Eq.(15) for 



on Eq.( 16 ) were presented to confirm the viability of our 
approach. Their straight forward extension to three di- 
mension thus promises to provide us with new insight on 
the fate of the heavy quarkonia, such as J/'^ and T in 
the hot QCD plasma as a function of real-time. 

In closing we list several open points of interest as well 
as possible further improvements: 

(i) The color charge is completely ignored in our de- 
scription. This is understood as a result of trac- 
ing out the color degrees of freedom in the master 



the density matrix Pqq (r,r',t) in the relative co- 
ordinates. However, it would be more appropriate 
to give a description of a wave function with color 
degrees of freedom and derive a more direct connec- 
tion to QCD. 

(ii) In addition to quantum decoherence, which we 
study in this paper and whose classical counter- 
part is the noise term in the classical Langevin the- 
ory, the quantum counterpart of friction is another 
important ingredient in the description of an open 
quantum system. This aspect will become essential 
in kinetic thermalization and the late time evolu- 
tion of heavy quark systems as well as resolve the 
question regarding the monotonous energy increase 
observed in our numerical simulation. 

(iii) A first principles definition of r(X, X') is desirable. 
Physical intuition have lead us to relate the thermal 
wavelength /th and the correlation length of the fluc- 
tuation. Since it gives the length scale of random 
fluctuation, the entropy change induced by insert- 
ing a heavy quarkonium might thus be related to 
T{X,X'). It would also be interesting to discuss 
thermodynamic quantities, such as free energy, in 
terms of our description. 

(iv) Application to relativistic heavy ion collision is in- 
teresting and possible in principle, but it would 
require some phenomenological refinements, e.g. 
how to combine non-relativistic quantum mechanics 
with the relativistic hydrodynamic expansion, how 
to describe hadronization of heavy quarkonium or 
heavy mesons by using the wave function at the 
freezeout, and so on. 
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